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Impact of neutron star oscillations on the accelerating 
electric field in the polar cap of pulsar 

or could we see oscillations of the neutron star after the glitch in pulsar? 
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Abstract Pulsar "standard model", that considers a 
pulsar as a rotating magnetized conducting sphere sur- 
rounded by plasma, is generalized to the case of oscil- 
lating star. We developed an algorithm for calculation 
of the Goldreich- Julian charge density for this case. We 
consider distortion of the accelerating zone in the polar 
cap of pulsar by neutron star oscillations. It is shown that 
for oscillation modes with high harmonic numbers {I, m) 
changes in the Goldreich- Julian charge density caused by 
pulsations of neutron star could lead to significant alter- 
ing of an accelerating electric field in the polar cap of pul- 
sar. In the moderately optimistic scenario, that assumes 
excitation of the neutron star oscillations by glitches, it 
could be possible to detect altering of the pulsar radioe- 
mission due to modulation of the accelerating field. 

Keywords starsmeutron, oscillations, magnetic fields • 
pulsars: general 



1 Introduction 

Neutron stars (NS) are probably the most dense objects 
in the Universe. There are extreme physical conditions 
inside NS, i.e. the magnetic field is close to the quantum 
limit, the pressure is of the order of the nuclear one and 
the typical radius of a NS is only about 2-3 times larger 
than its gravitational radius. Knowledge of properties of 
matter under such extreme conditions if very important 
for fundamental physics. It is impossible to reconstruct 
such physical circumstances in terrestrial laboratories, 
therefore, study of the NS's internal structure would give 
an unique opportunity for experimental verification of 
several fundamental physical theories. 

To study interiors of a celestial body one have to 
perform some kind of seismological study, by compar- 
ing observed frequencies of eigenmodes with frequen- 
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cies inferred from theoretical considerations. Eigenfre- 
quencies of oscillations in the crust of NS as well as in 
its interiors were calcti l ated in several p apers (see e.g. 
iMcDermott et al ] ll988t lGhugunovll200 fi( and references 
there). However, for seismological study there should ex- 
ist both i) a mechanism for excitation of oscillations, ii) a 
mechanism modulating radiation of the celestial object. 

There are two types of known NSs: member of binary 
systems and isolated ones. The former radiate due to ac- 
cretion of the matter from the companion. For these stars 
there are many possibilities to excite oscillations, for ex- 
ample by instabilities in the accretion fiow. However, in 
this case it would be difficult to distinguish whether a 
particular feature in the power spectrum of the object 
is due to oscillations of the NS or it is caused by some 
processes in the accretion disc/column. Because of this 
ambiguity we think that the study of isolated NS should 
be more promising in regard of the seismology. 

The vast majority of known isolated neutron stars are 
radiopulsars. The glitch (sudden change of the rotational 
period) is probably the only possible mechanism for ex- 
citation of oscillations for isolated pulsars. Radiation of 
radiopulsars is produced mostly in the magnetosphere. In 
order to judge whether oscillations of the NS could pro- 
duce detectable changes in pulsar radiation, the impact 
of the oscillations on the magnetosphere must be consid- 
ered. There is a widely accepted model of radiopulsar as a 
highly magnetized NS sur rounded by non-neutral plasma 
l|Goldreich fc Juhanl IT969.) . Although, there is still no 
self-consistent theory of radiopulsars, there is a general 
agreement regarding basic picture for the processes in 
the magnetosphere. Oscillations of the star can generate 
electric field as it happens in the case of rotation. Gen- 
eralization of the formalism developed for rotating NS 
to the case of oscillating star should help to obtain the 
desired information. 

The first attempt to generalize Goldreich- Julian (GJ) 
formalism to the case of oscillating NS was made in 
iTimokhin et al. I (|200d). It was developed a general al- 
gorithm for calculation of the GJ charge density in the 
near zone of an oscillating NS. Using this algorithm GJ 
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charge density and electromagnetic energy losses were 
calculated for the case of toroidal oscillations of the NS. 
Here we apply this formalism to the case of spheroidal os- 
cillation modes, representing wide class of stellar modes 
(r-,g-,p- modes). We consider also impact of stellar oscil- 
lations on the acceleration mechanism in the polar cap of 
pulsar and discuss the possibility of observation on the 
NS oscillations. 



2 Main Results 

2.1 General formalism 

Let us start by considering the case of a non-rotating os- 
cillating NS. Motion of the conducting NS surface in the 
strong magnetic field of the star generates electric field 
as in the case of rotation. Only oscillation modes with 
non-vanishing velocity V""" at the surface will disturb 
the magnetosphere. For the same reason as it is in the 
pulsar "standard model", the electric field in the mag- 
netosphere of an oscillating star should be perpendicular 
to the magnetic field. Otherwise charged particles will 
be accelerated by a longitudinal (parallel to B) electric 
field and their radiation will give rise to electron-positron 
cascades producing e nough particles to screen the accel- 
erating electric field l|Sturrocklll97l|) . As in the case of 
rotating stars we will define the Goldreich- Julian electric 
field Egj as the field which is everywhere perpendicular 
to the magnetic field of the star, Eqj -L B, and the GJ 
charge density as a charge density, which supports this 
field 
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For simplicity we consider only a zone near the NS, 
at the distances r <ti i-Kc/uj, where oj is the frequency of 
NS oscillations. For rnany g lobal oscillation modes (see 
e.g. iMcDermott et aP igSS*) the polar cap accelerating 
zone is well withing this distance, therefore, we can study 
the changes of the accelerating electric field in the polar 
cap caused by oscillation. In the near zone all physical 
quantities change harmonically with time, i.e. the time 
dependence enters only through the term e^"^*. 

We make an additional assumption, that changes of 
the magnetic field induced by currents in the NS crust 
are much larger than the distortion caused by currents 
fiowing in the magnetosphere. This assumption is con- 
sidered as a first order approximation according to the 
small parameter (f/i?Ns)j where ^ is the amplitude of os- 
cillation and i?Ns~ the NS radius. In other words, outside 
of the NS 



V X B = 



(2) 



in the first order in (^/i?Ns)- This assumption can be 
rewritten in terms of a condition on the current density 
in the magnetosphere as 
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Condition ||3J| implies that the current density in the near 
zone of the magnetosphere is less that the GJ current 
density connected with oscillations multiplied by a large 
factor c/{ujr). So, if the current density in the magneto- 
sphere is of the order of the GJ current density, assump- 
tion JSJ is valid. Under these assumption it is possible 
to solve the problem analytically in general case, i.e. to 
develop an algorithm for finding an analytical solution 
for the GJ charge density for arbitrary configuration of 
the magnetic field and arbitrary velocity field on the NS 
surface. 

Under assumption ^ the magnetic field can be ex- 
pressed through a scalar function P as 



B = V X V X (Per) 



(4) 



The GJ electric field depends also on a scalar potential 
'Fqj through the relation 



Egj = --V X (dtPe^) - Vifc,, • 
c 

The GJ charge density is then expressed as 

Pg.i = -l^^'^Gj • 
47r 

An equation for S'g, is 



(5) 
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{drd^P dgdtP - drdeP O^OtP) - , (7) 



where Aq is an angular part of the Laplace operator. 
This is the first order linear partial differential equation 
for the GJ electric potential iJ'gj- As the equation for 
iJ'gj is linear, each oscillation mode can be treated sepa- 
rately. This equation is valid for arbitrary configuration 
of the magnetic field and for any amplitude of the sur- 
face oscillation, provided that condition ^ is satisfied. 
Dependence on oscillation mode appears in the bound- 
ary conditions and also through the time derivative dtP- 
For different oscillation modes both the equation and 
boundary conditions for tf^cj are different. Derivation of 
eq. (0, boundary condition for functions !Z'gj as well as 
more de tailed discussion of use d approximations can be 
found in iTimokhin et al.l Il2n0f]t) . 



2.2 Goldreich Julian charge density 



(3) 



In lTimokhin et all (|200(tl equation was solved for the 
case of small-amplitude toroidal oscillations and dipolar 
configuration of unperturbed magnetic field. Solutions 
had been obtained with a code written in computer alge- 
bra language MATHEMATICA. Now we have developed 
a new version of this code, which allows to obtain analyt- 
ical solutions of equation Q for a more complicated case 
of spheroidal modes. Any vector field on a sphere can be 
represented as a composition of toroidal (V ■ V""" = 0) 
and spheroidal (V x Y'^" — 0) vector fields l|Unno et al.l 
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Fig. 1 The shape of the star surface during oscillations left for spheroidal mode (7, 3), right for spheroidal mode (7, 2). The 
shape changes harmonically with time. 




Fig. 2 Electric potential •f'c.j along a dipolar magnetic field line as a function of the polar angle 6 is shown for 5 field 
lines with azimuthal angle (j) ~ 0- l^ft for spheroidal mode (7,3), right for spheroidal mode (7,2). The quantity changes 
harmonically with time. 



So, now we are able to calculate GJ electric field 
and charge density for arbitrary oscillations of a NS with 
dipole magnetic field. 

Similar to the case of toroidal oscillations the small 
current approximation turned to be valid for a half of 
all oscillation modes. For oscillation modes with veloc- 
ity field, which is symmetric relative to the equatorial 
plane (see an example of such mode in Fig. ^ (left)), so- 
lution of eq. (UI) is smooth everywhere (see Fig. |3 (left)). 
For modes with antisymmetrical velocity field (an ex- 
ample of such mode is shown in Fig. ^ (right)), iJ'gj is 



discontinuous at the equatorial plane (see Fig|21 (right)). 
There is the following reason for such behavior. Dipolar 
magnetic field is antisymmetric relative to the equato- 
rial plane. Antisymmetric motion of the field line foot- 
points give rise to a twisted configuration of the magnetic 
field, which cannot be curl-free. So, for such modes a 
strong current will flow along closed magnetic field lines, 
j 3> Jg-i- However, there is no physical reason why a 
smooth solution for GJ electric field can not exists also 
for such modes. An argument supporting this hypothesis 
is a solution for twisted force-free magnetic field found 
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Fig. 3 Charge density pcj near tlie NS for azimuthal angle 
and negative ones in felue. left: pcj for the spheroidal mode 
changes harmonically with time. NB: at the equatorial plane 




(p = 0. Positive values of the charge density are shown in red 
(7,3); right: pcj for the spheroidal mode (7,2). The quantity 
{9 = 7r/2) Pg'j is infinite. 



Fig. 4 Charge density pGj{r = RNs,6,(t> = 0) on the NS surface is shown in a polar coordinate system ( |pg^(-Rns, 6', 0) | ,9). 
Positive values of charge density are shown by the solid line and negative ones by the dashed line, left: pcj for the spheroidal 
mode (7,3); right: pcj for the spheroidal mode (7,2). The quantity changes harmonically with time. Circles correspond to 
the values of poj = (0.1, 0.2, 0.3, 0.4, 0.5) in normalized units. NB: on the equatorial plane {6 = it/2) pcj is infinite. 



bv lWolfsoiJ l[l99?l. In his solution, which corresponds to 
the toroidal mode (2,0), the configuration of force- free 
twisted magnetic field is supported by a strong current 
flowing along magnetic field lines. 

For the modes with smooth solutions our approxi- 
mation should be valid. On the other hand, a strong 
electric current will flow only along closed magnetic fleld 
lines. The current density along open magnetic field lines 
should be close to Jg.t and condition for the small cur- 
rent approximation (eq. Q), will be satisfied in the open 
field line domain. Then, Wq^ in the polar cap could be ob- 
tained by solving eq. (0. ^gj in the polar cap will differ 
from the solutions obtained here, because the boundary 
conditions should be set at the polar cap boundaries and 
not on the whole surface on the NS. The main proper- 
ties of regarding its qualitative dependence on the 



coordinates, used in our discussion, would be, however, 
similar to the properties obtained from our solutions. 

As expected, the GJ charge density distribution fol- 
lows the distribution of the velocity field (see e.g. Figs.|3IElEl). 
With increasing of the harmonic numbers, pcj tsills more 
rapidly with the distance, what is also expected for mul- 
tipolar solutions. A remarkable property of GJ charge 
density distribution near oscillating star is that the local 
maxima of pcj increases with increasing of both I and m 
(see Figs. El IZ|)- The reasons for this are as follows. The 
electric field induced by oscillations is of the order 

Eg., ~ —B . (8) 
c 

The charge density supporting this electric field is of the 
order of E/Ax, where Ax is a characteristic distance of 
electric field variation. For a mode with harmonic num- 
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Fig. 5 View of the polar cap of pulsar from the top. 
Goldreich- Julian charge density for oscillation mode (44, 4) 
is shown by the color map (positive values in red, negative - 
in blue). The quantity changes harmonically with time. The 
polar cap boundary for a pulsar with period 3 ms is shown 
by the dashed line. 



ber I this size is of the order of R^s/i- Hence, for pc, we 
have an estimate 



I 



E 



B 



4nR^ 



(9) 



So, for the same amplitude of velocity of oscillation the 
amplitude of variation of the GJ charge density is larger 
for higher harmonics. Equation (jS)) is in a good agree- 
ment with the exact results shown in Figs. [7| 



2.3 Particle acceleration in the polar cap of pulsar 

The GJ charge density induced by NS oscillations influ- 
ences particle acceleration mechanism in the polar cap 
of pulsar. As we will show, oscillations will have the 
strongest impact on the accelerating electric field in the 
polar cap of pulsar for models with free particle escape 
from the NS surface ijScharlemann et al.ll978HMuslimov ^ 
The accelerating electric field in pulsars arises due 
to deviation of the charge density of the plasma from the 
local GJ charge density. For pulsar models with Space 
Charge Limited Flow (SCLF) the charge density of the 
flow p at the NS surface is equal to the local value of the 
GJ charge density, p(i?Ns) = Pgj(-Rns)- Magnetic field 
lines diverge and the charge density of the flow decreases 
with increasing of the distance from the star. However, 
the local GJ charge density decreases in a different way 
and at some distance r from the NS Pg.j(^) 7^ Pi''')- The 
discrepancy between these charge densities gives rise to 



a longitudinal electric field (see lScharlemann et alJUQTSt 
iMusHmov fc TsvganlllQQl . 

The GJ charge density for oscillation modes with 
large /, m falls very rapidly with the distance. Hence, 
the charge density of a charge-separated flow for oscil- 
lating NS will exceed the local GJ charge density at some 
distance from the star. This produces a decelerating elec- 
tric field. Therefore, in the case of rotating and oscillat- 
ing NS, if pj^fj < pJJj, the effective accelerating electric 
field will be reduced periodically due to superposition of 
accelerating and decelerating electric fields. 

The oscillational GJ charge density p"^^ for modes 
with large I, m decreases practically to zero already at 
2-3 NS radii and the whole oscillational GJ charge den- 
sity contributes to the decelerating electric field. While 
in the case of rotation only ~ 15% of the GJ charge den- 
sity contributes t o the accelerating electric fie ld for SCLF 
in the polar cap l|MusHmov &: TsvgaiJll992() . The most 
important factor increasing modification of the acceler- 
ating electric field is that the amplitude of j increases 
with increasing of the harmonic numbers of the mode. In 
Fig. Owe show the difference between the charge density 
of SCLF and the local GJ charge density for dipolar mag- 
netic field as a function of the distance for rotating and 
oscillating stars. In order to demonstrate the importance 
of the discussed effects we consider a rather grotesque 
case when the linear velocity of rotation is equal to the 
maximum velocity of oscillations. It is evident from this 
plot, that for large enough I and m the decelerating elec- 
tric field caused by stellar oscillations could be of com- 
parable strength with the accelerating electric field even 
if < V'°\ 

Let us estimate the harmonic number of the oscilla- 
tion mode where decelerating electric field would have 
a given impact on the accelerating electric field induced 
by NS rotation. The decelerating electric will be k times 
less that the rotational accelerating electric field. 



acc 



0.15pg 



if 



I - 0.15 k 



V" 



(10) 



(11) 



1E*9%:¥]J' is the linear velocity of NS rotation at the equa- 
tor. For example for I ~ 100, from this equation we can 
see, that for canceling of the accelerating electric field it 
is sufficient that the velocity amplitude is only ~ lO^'^ 
of rotational velocity at the NS equator. 

For pulsar operating in lRuderman fc Sutherland! ^197,'^^ 
regime the impact of stellar oscillations on the acceler- 
ating field in the polar cap will be reduced by the factor 
of ^ 10. In this model the accelerating electric field is 
generated in a vacuum gap, so the whole rotational GJ 
charge density contribute to the accelerating electric field 
and in this case k ^ Pgi/Pgj- 



6 



A. N. Timokhin 



2.5 





Fig. 6 Dependency of the Goldreich- Julian charge density on 
the harmonic number /. Distribution of poj on the NS sur- 
face in the polar cap of pulsar is shown for oscillation modes 
(from left to right) (28,4), (44,4), (64,4) in polar coordi- 
nates (|pgj(^?ns, ^, 0) I ,6). Positive values of charge density 
are shown by the solid line and negative ones by the dashed 
line. The red line shows the angle at which the last closed 
field line intersect the NS surface for a pulsar with period 3 
ms. Circular segments correspond to the values of pa J in nor- 
malized units shown at the left side of each plot. The same 
normalization for poj is used as in Fig. ^ 
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Fig. 7 Dependency of the Goldreich- Julian charge density 
on the harmonic number m. Distribution of pcj on the NS 
surface in the polar cap of pulsar is shown for spheroidal 
oscillation modes (from left to right) (64,4), (64,8), (64, 14) 

in polar coordinates (Jpg^(-Rns, ^, 0) | ,9). Meaning of the lines 
is the same as in Fig|bl 



3 Discussion 

We have shown, that oscillations of the NS can induce 
changes in the accelerating electric field, which are more 
stronger than a naive estimation {V°"'/c)B. Indeed, for 
high harmonics the induced electric field will be ^ I 
times stronger. In order to make definitive predictions 
about observable parameters is necessary to study the 




Fig. 8 Ap = p — Pgj - difference between the charge den- 
sity of a space charge limited flow and the local Goldreich- 
Julian charge density along magnetic field lines in the polar 
cap of pulsar. Ap is shown (in arbitrary units) by solid lines 
for spheroidal modes with different {l,m), from top to bot- 
tom: (64,14), (54,14), (64,2), (54,2). The same relation for an 
aligned rotator is shown by the dashed line. Negative values 
of Ap give rise to an accelerating electric field, negative ones 
- to a decelerating field. 



polar cap acceleration zone more detailed. An accelerat- 
ing electric potential and the height of the pair formation 
front should be calculated. 

However, we can do simple estimations using equa- 
tion (|ll|l . Let us estimate the harmonic number of the 
mode which can cancel the accelerating electric field, as- 
suming the mode is excited by a glitch. As we pointed in 
Sec. l2.1l onlv modes with non-zero amplitude on the NS 
surface can produce changes in the magnetosphere. The 
distribution of oscillational motion plays crucial role. If 
oscillations are trapped in the NS crust a rather small 
energy will be required to pump the oscillation ampli- 
tude to the level high enough for strong disturbance of 
the accelerating electric field. If a fraction e of the NS 
mass Mns is involved in the oscillations the amplitude of 
the oscillational velocity is of the order 



itch 



(13) 



where W°''° is the total energy of the mode. The energy 
transferred during the glitch of the amplitude Af2 is 

'27ry Ail 

where i is the fraction of the total momentum of inertia 
of the NS /ns coupled to the crust. Let us assume that 
some fraction ry of this energy goes into excitation of 
oscillations. Using eqs. ((T^ . ((T^ . ((TT|l we get conditions 
for harmonic number of modes which would periodically 
cancel the accelerating electric field, as 

-1/2 

(14) 



1 > 300 r;-;/' £1/2 



rAQ^ 

1^/6 

for SCLF model, and for lRuderman fc Sutherland! I^T^) 
model: 



1 > 2000 



1/2 ^1/2 --1/2 f^^\ 



-1/2 



(15) 
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Here ry% is measured in per cents and the relative magni- 
tude of the glitch {Ail/il)^ is normalized to 10~^. It is 
widely accepted, that the origin of pulsar glitches is an- 
gular momentum transfer from the NS core to the crust. 
In the frame of this model the fraction of the energy 
which can go into excitation of NS oscillation 77 is of 
the order of Af2/f2, i.e. it is very small, of the order of 
~ 10~^. We may speculate however, that excited oscil- 
lations are trapped in the NS crust, i.e. e is also very 
small. Such global oscillation modes (/ ~ several hun- 
dreds) could induce substantial changes in the accelerat- 
ing electric field. 

As we mentioned in Sec. 12. H all physical quantities in 
the solutions obtained here oscillate with the frequency 
of the star oscillations. The accelerating electric field 
close to the local geometrical maxima of the oscillational 
GJ charge density will be weakened periodically by the 
decelerating effect due to stellar pulsations. The field os- 
cillation will influence the particle distribution in the 
open field line zone of the pulsar magnetosphere and 
it should produce some observable effects. Depending 
on oscillation mode and position of the line of sight a 
complicated pattern will appear periodically in individ- 
ual pulse profiles. Although individual pulses are highly 
variable, the presence of periodical features should be 
possible to discover in the power spectra of pulsars, pro- 
vided the oscillations are excited to a high enough level 
and observations have been made with hight temporal 
resolution. If one observes some feature, which appears 
just after the glitch, then decreases and disappears after 
some time, and it never appears in the normal pulsar 
emission, then one can undoubtedly attribute this fea- 
ture to the NS oscillations. 
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